Hamiltonian derivation of the Charney-Hasegawa-Mima equation 
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The Charney-Hasegawa-Mima equation is an infinite-dimensional Hamiltonian system with dy- 
namics generated by a noncanonical Poisson bracket. Here a first principle Hamiltonian derivation 
of this system, beginning with the ion fluid dynamics and its known Hamiltonian form, is given. 

I. INTRODUCTION 

When dissipative terms are dropped, all of the important models of plasma physics are described by partial differen- 
tial equations that possess Hamiltonian form in terms of noncanonical Poisson brackets. For example, this is the case 
for ideal magnetohydrodynamics [110,0], the Vlasov-Maxwell equations 4, 5, 6|, and other systems (see Refs. 7, 8, 9| 
for review). Among these, there exist several reduced fluid models whose Hamiltonian structure has been derived a 
posteriori. These include the four- field model for tokamak dynamics of Hazeltine et _al. |10| : models for collisionless 
magnetic reconnection derived and investigated by Schep et al. Kuvshinov et al 
recent gyrofluid model of Waelbroeck et al 



12l |. and Tassi et al. 



13| ; and the 



14l | . The noncanonical Hamiltonian formulation has also been adopted to 
investigate the electron temperature gradient driven mode [15] and convective-cell formation in plasma fluid systems 
In addition to these fluid models, the Hamiltonian structure of kinetic and reduced kinetic equations has also 
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2l| for review). 



been highlighted, for example, in guiding-center theory and gyrokinetics (see Refs. 

This Hamiltonian form originates from the Hamiltonian and action principle forms of the basic electromagnetic 
interaction, i.e., the Hamiltonian form possessed by the equations that describe a system of charged particles coupled to 
Maxwell's equations (see, e.g., Ref. 9| for discussion). It is now well established that there exist numerous advantages 
of such a Hamiltonian formulation, among which are the identification of conserved quantities (that are important for 
the verification of numerical codes), the study of stability, the use of techniques for Hamiltonian systems like averaging 
and perturbation theory, etc. Here we perform a perturbative derivation within the noncanonical Hamiltonian context, 
which means the Poisson bracket as well as the Hamiltonian must be expanded. 

In a nutshell, a Hamiltonian system is a system whose dynamics of any observable F (depending on a finite or 



infinite number of variables) can be written using a Hamiltonian (scalar) function H and a Poisson bracket {•, •} as 

dF 

where the Poisson bracket satisfies the following properties: bilinearity, antisymmetry, Leibniz rule, and Jacobi 
identity. Given a reduced model whose dynamics is given by a partial differential equation, it is in general difficult to 
guess whether or not the model is a Hamiltonian system, and if it is, finding the Hamiltonian and the Poisson bracket 
may be similarly difficult. There are basically two methods for finding Hamiltonian structure: the first method is to 
use physical intuition to obtain the Hamiltonian (energy) and to construct a general class of antisymmetric operators 
which, when acting on the gradient of the Hamiltonian, produces the equations of motion. Then, the Jacobi identity is 
used to select from the class the desired operator that is the essence of the noncanonical Poisson bracket. This method 
has been used to obtain a large number of basic and approximate Poisson brackets for fiuid and plasma dynamics, 
examples being the reduced fiuid models cited above. The second method begins from a known or postulated action 
principle, in the latter case obtained by using physical intuition to obtain the 'energies' of the Lagrangian. Usually 
associated with the action principle is a canonical Hamiltonian description, which can be written by means of the 
chain rule in terms of physical variables of interest (e.g. Refs. [5|, Isi]) resulting in a noncanonical Poisson bracket. 

If one begins from some Hamiltonian parent model, some basic starting point in the derivation, and introduces crude 
approximations suggested, e.g., by physical considerations of some experimental set-up, then the Hamiltonian structure 
can be easily destroyed. The Hamiltonian form of the resulting system must therefore be verified, in particular, the 
Jacobi identity for the Poisson bracket. Given this verification, the reduced model is naturally equipped with a 
Hamiltonian structure since the Poisson bracket and the Hamiltonian function are provided by the derivation process 



(for an example of this derivation process, see Ref. 



m 



In this paper we consider the derivation of the Charney-Hasegawa-Mima (CHM) equation 
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both the dynamics of Rossby waves of geophysical fluid dynamics (see, e.g., Ref. 25|) and drift waves in inhomogeneous 



241 ■ which describes 



plasmas (see, e.g., Ref. 26|). We focus on the derivation in the plasma physics context but our analysis can be 
easily adapted to the geophysical context. In particular, we show how the Hamiltonian structure is preserved in 
the derivation of the CHM equation starting from a fluid parent model. In the present approach the Hamiltonian 
structure is provided by the derivation process, and the Jacobi identity need not be checked. 

We consider a plasma under the influence of a constant and uniform magnetic field B = Bi. The relevant dynamics 
occurs in the (two-dimensional) transverse plane whose coordinates in a given basis are denoted by x and y. Under 



some assumptions, the CHM equation gives the following evolution of the electrostatic potential (j){x,y,t) generated 
by the plasma: 

d 



dV 



A(/.)-[(/.,A0 + A] 



(1) 



where the bracket [■, ■] is given by 

df dg df dg . 

and A is any function of x and y (related to the equilibrium configuration). The infinite-dimensional phase space 
is composed of the variables 4){x,y) for any point {x,y) in the transverse plane. The space of observables, F, for 
this system is composed of functionals of (j>. It has been shown in Ref. ^] that this equation possesses an infinite- 
dimensional Hamiltonian structure where the Hamiltonian is 

and the noncanonical Poisson bracket is 



{F, G} = - J (fx{ct) - A0 - A) [(1 - A)-iF^, (1 - A)-iG^] 



(2) 



where F^f, denotes the functional derivative of the functional F with respect to the variable This Hamiltonian 
structure was found ad hoc in Ref. 27| by an educated guess in analogy with the vorticity equation for two-dimensional 

This analogy is rather straightforward if we consider the dynamics for the field 



a; 



incompressible fiow (see e.g. Ref. 
q — A(/) — + A which is given by the Hamiltonian 



and the Lic-Poisson bracket 



H = ]^j d^x{q-X}(l-A)-\q-X), 



{F,G}= J d^xq[Fg,Gg], 



J| and a quite general class of systems 



281. In what 



which is of the same form as that for the Vlasov-Poisson system 
follows, we start by considering a Hamiltonian formulation for the fiuid equations for the ions (in Sec.[Tl|) and derive 
the above Hamiltonian and Poisson bracket from this formulation (in Sec. IIIip . 
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II. ION FLUID EQUATIONS AS A HAMILTONIAN SYSTEM 



We start the derivation of the CHM equation from two dynamical equations: one describing the transverse dynamics 
of the ion velocity field v(a;, y, t) and the other describing the dynamics of the ion density field n{x, y, t): 

M (v + (v • V)v) = -eVcj) + ev X B, (3) 
7T, = -V • (nv), (4) 

where the dot indicates the partial derivative with respect to time t. The electrostatic potential (j) is obtained from 
the dynamics of the electrons: by neglecting their inertia, the electron density obeys the Boltzmann law 



rie — HQ exp {e4>/T) . 



(5) 



where T is the electron temperature and no = nQ{x, y) is the electron density at equilibrium. From the quasi-neutrality 
condition, we obtain that n — n,,. The total energy of the ions, given by the sum of their kinetic energy plus the 
potential energy provided by the electric field, is a conserved quantity that is also a good candidate for the Hamiltonian 
of the system of Eqs. ([3]|4]). This Hamiltonian is written as 

(6) 



H{n,w) 

The dynamics is determined by the Poisson bracket 
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{F.G} 



d^x 



■ VGn — Gv • V-Fn 



V X 



- + — z • Fv 



xGv) 



(7) 



where ujc — eB/M. The bracket of ([7]) is identical to a portion of that of Ref. 1| with the inclusion of an additional 
'vorticity' term, LUcz/n; consequently, it is known to satisfy the Jacobi identity. It is easy to verify that the ion 
momentum equation is obtained from the bracket of the velocity field with the Hamiltonian ([6]): 

V = {v,H} = -(v • V)v - — Vln f — I +u}cv x z, 

M ynaj 

and, similarly, the ion continuity equation is given by 

h = {n, H} = —V ■ (nv). 

III. CHARNEY-HASEGAWA-MIMA EQUATION 



Without loss of generality we write the vector field v{x, y, t) in terms of two scalar fields </) and T as 

v = z X V0 + VT, 



(8) 



where one function is related to V • v and the other to V x v by the relations: = z • V x v and AT = V ■ v. In 
fact, we find it more convenient to consider a related change of variables (n, v) (—> (n, q, D) defined by 

n — n, 

Z • V X V + LJc 



q = 



n 



£> = Vv. 

The above equations are incomplete because they do not possess a unique inverse. However, a unique inverse is defined 
by the following: 

n — n, 

V = z X VA~^{qh - uic) + VA^^D, 



where 



A-^F = ^— /d2a;'ln||x-x'||i^(x'). 
27r J 



In terms of the new variables (h^q, D), the Hamiltonian ([6]) becomes 

|VA-i(gn-u;e)l' 



H{n,q,D) = / d^x 



[A-\qh-uj,),A-^D] 



\yA-^D\ 



since |z x V/p = iV/j'^ for any function f oi x and y, and the bracket ([7]) becomes 

{F, G]^ - ! (fx {-VFd ■ VGs + VGd ■ VFs + ^VFd ■ Vg - ^VGd ■ 
J \ n n 

qIFd^Gd] 



Fq Gq 



(9) 



It should be noted that Casimir invariants of such a bracket, which are the functionals that Poisson commute with 
all the other functionals ({C, G} — Q for all functionals G), are given by 



C = / (fxnT{q), 



where J- is any function of q. 

We first assume that the variables evolve slowly with time, which is equivalent to adding a factor of 1/e in front of 



the Hamiltonian, 



H{n,q,D) = - cfx 



+- 



— n fin I — 
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then we introduce an e-ordering for the dynamical variables. The hypothesis is that the system of interest is near an 
equilibrium state whose spatial variations are of order e: 

n(x, t) = no (ex) + eui (x, t) , 
v(x,t) = evi(x,i), 

which translates into an assumption on the new variables (n, q, D) and, in particular, on the definition of new dynamical 
variables (ni,gi,£>i), 

n = no + efii, 
q = qo + eqi, 
D = eDi, 

where qo = ojc/fio and ng = no (0,0) are constant (the spatial variations of no are included in ni). Notice that the 
potential energy can be rewritten as 

no 



n In — , 
no 



with the following expansion: 



fi(ln( ^) -l) =-no-noln^ + e2^-eniln^ + 0(e3). 
no / / no 2no no 



The term — en-i ln(no/no) is of order e^, due to the spatial variations of no, which can be seen by writing no = fio + edno: 

_ no 2"^l<5^0 , 3^ 

— eni In — = — e h 0(e ). 

no no 

Next, we expand the Hamiltonian and the Poisson bracket: the Hamiltonian is 



■I- 



H = e (fx ho 



|VA-i(gino + gorii)P ^ [VA-^Dip ^ T h\-2h^8no 



2M h?, 



since / d^a;no[A ^(gino + g'o'^i), A ^Di] = 0, and the Poisson bracket is 
{F,G} = d^xiVFo, ■ VGn, - VFn, ■ VGdJ 
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Thus the dynamics emerges to leading order at (which gives the dynamics on a time-scale of order e) and to next 
order at e° (whose influence happens on a time-scale of order one). 

First we study the dynamics given by the leading order. It should be noticed that qi is constant, since the leading 
order Poisson bracket does not contain any functional derivatives with respect to qi , and that 

= ^AHd, + 0(e) = -"^Di + [ni, A-i(gino + ^0^1)] - V • {hiVA-^Di) + 0(e), 

D, = -\AHf,, + -V- (^S7q,)+l[Hn„qi]+0{e). 
e V no / e 

If we impose the following constraints on the initial conditions: 

AFft, = 0, 
AHd, = 0, 

then these constraints are preserved by the leading order flow. These constraints are equivalent to the following: 

£'1=0, 

T Aril . , - - X „ 

■Tj— qono[qoni + qmo) = 0. 

M no 

Note that, from expanding n(x, t) — no(ex) + eni(x, t) about e = 0, it follows that Stiq is a linar function of x and, as 
a consequence, it does not appear in the equations for the constraints. A generalization to the case of non-harmonic 
6no is however possible. 

Even if wc neglect higher order terms (e^ in the Hamiltonian and in the Poisson bracket), these constraints are 
not preserved by the Poisson bracket. Therefore, these quantities are approximately conserved on a time scale of 
order e. Next, we approximate the dynamics on a time-scale of order 1 by inserting the constraints on Di and ni into 
the second order Poisson bracket. By dropping all dependence on Di and rii, the dynamics is thus equivalently given 
by the Hamiltonian 

"0 / / ~ , T h\-2ni6nQ 



i/i = y" <fx^ ^-(go^i +g'ino)A ^{qohi + qiho) 



M hi 



where fii is a function of qi given by 



no 

m = 



and the Poisson bracket 



/ T 










j (fxqi 








. no 


no _ 



{F,G}^= d'xq, , (10) 



8 

which satisfies the properties of a Poisson bracket - in particular, the Jacobi identity. Using this condition on ni, the 
Hamiltonian Hi can be rewritten as 



TflQ 

til = „ , , n fax 



(11) 



where A contains the spatial variations of the equilibrium density as follows: 

A = — — (5no(x). 

no 

Using the symmetry of the operator (1 — {T/Muj'^)A)~^, the Hamiltonian ITTI) can be rewritten as 

Trio f ,2 , .^ f. T 



Up to some constants, the Poisson bracket (fTUl) and the Hamiltonian (fT^ are indeed the same as those presented in 
Ref. ^2?!. Thus we have provided a derivation process that leads to dynamics, on time-scales of order 1, that is still 
generated by a Hamiltonian and a Poisson bracket. 

IV. CONCLUSION 

An important issue in the derivation of reduced models for plasma physics is avoiding the introduction of fake 
dissipative terms, which may result from uncontrolled approximations and truncations in the derivation process. In 
particular, if the parent model has a Hamiltonian structure, we argue that the final reduced model should also have 
a Hamiltonian structure. 

In this paper we examined, in this spirit, the case of the CHM equation. In particular we showed how the fun- 
damental elements, i.e., the Hamiltonian functional and the Poisson bracket, of the Hamiltonian formulation of the 
CHM equation, emerge from the Hamiltonian structure of a parent model, which is the starting point of the deriva- 
tion commonly adopted in the plasma physics literature. The appearance of the Hamiltonian and the bracket of the 
CHM equation in the derivation process was seen to be facilitated by adopting the new set of variables {q,h,D). In 
terms of these variables, the part of the bracket of the parent model that becomes the CHM bracket can be easily 
identified. Indeed, what our paper shows is how the ordering adopted in the derivation is able to reduce the bracket 
of the parent model to the CHM bracket, without compromising the fundamental properties of a Poisson bracket, 
such as for instance the Jacobi identity. A further new element of our analysis is the way the plasma compressibility 
is treated. Without invoking the drift approximation and the polarization drift, the divergence-free condition on the 



9 

plasma velocity appears as a solution for the variable Di on a time scale of order e. Such a solution is used in order 
to approximate the dynamics of order 1, assuming for such dynamics that an equilibrium solution exists. A similar 
argument is used for the dynamics of ni, which at the lowest order is constrained to be a function of qi or, more 
precisely, to be proportional to the plasma stream function. 

We believe that the method adopted in this paper is a framework for deriving the Hamiltonian structure in other 
reduced models of plasma physics, and for deriving new models while avoiding the risk of introducing fake dissipative 
terms. 
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